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ABSTRACT 

The TRACER cosmic-ray detector, first flown on long-duration balloon (LDB) in 2003 for 
observations of the major primary cosmic-ray nuclei from oxygen {Z — 8) to iron (Z = 26), 
has been upgraded to also measure the energies of the lighter nuclei, including the secondary 
species boron (Z = 5). The instrument was used in another LDB flight in 2006. The properties 
and performance of the modified detector system are described, and the analysis of the data 
from the 2006 flight is discussed. The energy spectra of the primary nuclei carbon (Z = 6), 
oxygen, and iron over the range from 1 GeV amu~^ to 2 TeV amu~^ are reported. The data for 
oxygen and iron are found to be in good agreement with the results of the previous TRACER 
flight. The measurement of the energy spectrum of boron also extends into the TeV amu~^ 
region. The relative abundances of the primary nuclei, such as carbon, oxygen, and iron, above 
^ 10 GeV amu^^ are independent of energy, while the boron abundance, i. e. the B/C abundance 
ratio, decreases with energy as expected. However, there is an indication that the previously 
reported dependence of the B/C ratio does not continue to the highest energies. 

Subject headings: Astroparticle physics — Cosmic rays — Balloons 



1. Introduction 

Direct measurements of the energy spectra 
of cosmic ray nuclei that extend into the high- 
energy region, i.e. energies around and above 
10^^ eV amu~^, require the exposure of large de- 
tectors over extended observation times. Only 
one such measurement has thus far been real- 
ized in space flight, the Cosmic-Ray Nuclei de- 
tect or (CRN) that was flow n on the Space Shut- 
tle (|L'Heureux et al. Il990l) . Long-duration bal- 
loon (LDB) flights permit exposures of the order 
of weeks, and several LDB instruments have been 
developed: TRACER (Transition Rad iation Ar- 
ray for Cosmic Energetic Radiation) (lAve et al 



20111 ). which is the subject of this paper, and two 
instruments that use short c alorimeters for the 
energy measur ement, ATIC ([Ganel et al.l 120051 ) 



and CREAM (AhnetaL 20071). TRACER was 
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designed and constructed at the University of 
Chicago, utilizing the heritage of CRN. In contrast 
to the calorimetric measurements, a nuclear inter- 
action is not required to occur, because charge 
and energy of cosmic rays are determined exclu- 
sively via electro-magnetic interactions. This al- 
lows TRACER to employ light-weight materials 
that consequently allow the instrument to be very 
large. It has a geometric factor of ~5 m^ sr (de- 
tector area 2x2 m^) and is currently the largest 
balloon-borne cosmic-ray detector. 

TRACER had a two-week circum-polar LDB 
flight in Antarct ica in 2003 (LDB l). The results 
from this flight ( Ave et al. 2008() provide a com- 
prehensive set of measurements of the individual 
energy spectra of the major primary cosmic-ray 
nuclei from oxygen {Z = 8) to iron {Z = 26), 
covering the energy range from ^ 1 GeV amu^^ 
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to several TeV amu"^. In total energy, the mea- 
surement reached energies well above 10^^ eV for 
the more abundant species. The results from 
this flight indicate that power-law fits in en- 
ergy to the measured spectra can be made above 
20 GeV amu^^ with the same power-law expo- 
nent, a = 2.67 ± 0.05, for all elemental species. 
There seems to be no noticeable dependence of 
this exponent on nuclear charge Z. This result 
appears to be surprising, and in fact, an anal- 
ysis of the TRACER results in the context of 
simpl e propagation a nd source models of cosmic 
rays (|Ave et al.ll2"ooi ) indicates that the observed 
spectra might not be exact power-laws but exhibit 
a slight concave curvature due to the competition 
of interaction losses and energy- dependent diffu- 
sive escape during the propagation through the 
Galaxy. If this curvature existed, it could be 
obscured by the statistical uncertainties in the 
measurement. 

Within measurement uncertainties, the TRACER 
data appear to be consistent with the results re- 
ported from ATIC and CREAM. In particular, 
the CREAM group also has reported a com- 
mon power law ex ponent of 2.66 ± 0.04 for the 
measured spectra ( Ahn et al. 2009f) . More re- 



cently, the CREAM group has claimed a "dis- 
crepant" upturn in the measured spectra at 
- 200 GeV amu-i (|Ahn et all boiof l but again, 
the statistical significance is quite limited. 

The cosmic-ray intensities measured near the 
Earth are significantly affected by Galactic prop- 
agation and do not directly reflect the charac- 
teristics of the cosmic-ray sources. At low en- 
ergies (< I GeV amu^^), the time for diffusive 
propagation through the Galaxy is ~ 1.5 • 10'' 
years, a t an interstellar path length of a few 
g cm ~^ ( Garcia- Munoz et al.lll977l : lYanasak et al. 
20011 ). It has been known for a long time that 
the path length decreases at relativistic ener- 
gies, appa rently as a power law proportional 
to E-°-^ (lJuliusson et al.l Il972t lYanasak etlir 



20011 lEngelmann et al.lllQQotlSwordv et al.lll99ol) 

Whether this strong decrease pe rsists to the high- 
est energies, is an open question ( Ave et al1[2009 '). 
While the particles diffuse through interstellar 
space, they can also be lost by spallation reac- 
tions with the components of the interstellar gas. 
The spallation cross sections depend on the mass 
number A of the projectile approximately like 



A^^^, hence heavy nuclei are more strongly af- 
fected than light ones. In order to determine the 
characteristics of the cosmic-ray sources, these 
two propagation effects, one depending on energy 
and the other on nuclear mass number, must be 
understood to deconvolve the measured spectra. 

The propagation path length of cosmic rays can 
naturally be determined by measuring the inten- 
sities and energy spectra of the secondary nuclei, 
i.e. those elements that are produced exclusively 
by interstellar spallation. This group includes 
the light elements lithium, beryllium, and boron, 
where boron is the most studied element. At 
relativistic energies, quite accurate measu rements 
exist from HEAO (Engelmann et al.lll990l) , but do 
not extend beyond ~ 30 GeV amu~^. At higher 
energies, up to a few hundred GeV amu ~^ the 



boron intensity has been measured by CRN (jSwordv et al 
19901) and results on the relative abundance of 



boron have been reported by CREAM ()Ahn et al 



2008), but the data are severely limited by statis- 
tical uncertainties at the highest energies. 

The TRACER instrument could not include an 
energy measurement of the light nuclei below oxy- 
gen in LDBl because of limitations in the dynamic 
range of its detector electronics. A major upgrade 
of TRACER was performed after that flight. The 
instrument is now sensitive to the entire charge 
range A < Z < 28 and was launched on an LDB 
flight from Kiruna, Sweden in 2006 (LDB2). In 
the following, we will briefly describe the mod- 
ified detector system, discuss the data analysis, 
and present the results from the this flight, as well 
as the combined data set from LDBl and LDB2. 

2. The Measurement 

2.1. The Detector and Measurement Prin- 
ciple 

The TRACER instrument represents a combi- 
nation of electronic detectors, comprising plastic 
scintillator sheets, acrylic Cerenkov counters, and 
arrays of single- wire proportional tubes. The over- 
all dimensions of the detector system are quite 
large, leading to an aperture of 4.73 m^ sr. The in- 
strument and the measurement process for LDBl 
have b een described in some detail bv lAve et al 



(|2008[ ). and a broader, technical description which 
also includes the modification s for LDB2, is be- 
ing published ( Ave et al]l2011 ). Therefore, we will 



2 



Cerenkov 
Scintillator 
dE/dX Array 




Fig. 1. — Schematic view of the detector elements 
of TRACER. 



here only briefly review the measurement process 
and emphasize some aspects which are of specific 
importance for the 2006 measurement. 

A schematic drawing of TRACER is shown in 
Figure [TJ There are two identical pairs of scin- 
tillator and Cerenkov counters on top and bot- 
tom of the instrument, respectively, with each 
counter 200 cm x 200 cm in area. Each counter 
is viewed by 24 photomultiplier tubes (PMTs) via 
wavelength-shifter bars. The scintillators provide 
the coincidence trigger for the measurement. 

The space between these counter-pairs is oc- 
cupied by arrays of gas-filled, single-wire propor- 
tional tubes and layers of plastic-fiber radiators. 
There is a total of 1584 tubes, each 2 cm in 
diameter and 200 cm in length. They are ar- 
ranged in double-layers in densest packing and 
are alternately oriented in two orthogonal direc- 
tions. Thus, there is a total of 8 double layers of 
2 X 99 tubes, each. The four upper double lay- 
ers measure the ionization energy loss ("dE/dx"). 
Each of the four remaining layers is located be- 
low a radiator layer and measures ionization loss 
with transition radiation superimposed ("dE/dx 
-I- TR"). The arrangement and packing of the 
fiber radiators is i dentical to that used for the 
'Heureux et al. 1990l ). and for 



CRN instrument 



LDBl (|Ave et al.l 20081 ). A cosmic-ray particle 
passing through the instrument will traverse 16 
individual tubes and generate 16 independent sig- 
nals along its trajectory. The gas used in LDB2 
is a 95%:5% Xe:CH4 mixture at a pressure of one 
atmosphere. 

For each cosmic-ray nucleus that traverses the 
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Fig. 2. — Response functions, in units of sig- 
nal at minimum ionizing energy (MIP) normalized 
by Z^, of sub-detectors used for energy measure- 
ments. The dashed lines indicate the signals at 
minimum ionizing energy and TR onset. 



TRACER instrument, the nuclear charge and 
the energy E (or the Lorcntz- factor 7 — E/mc^), 
as well as the geometric trajectory through the 
instrument are determined. 

The charge is obtained from the signals of the 
Cerenkov-scintillator combination, utilizing the 
fact that electromagnetic interactions scale with 
Z^. However, the energy dependence of these 
two signals is different: the scintillation signal de- 
creases with /3~^, and the Cerenkov signal rises 



as (1 — 1/n P ) above a threshold /3 



(with 



refractive index n=1.49, and /? = v/c). With 
/3 ~ 1, the scintillation signal becomes indepen- 
dent of energy at a "minimum ionization level." 
Thus, if the two signals are combined, the depen- 
dence on energy (or /3) can be eliminated, and 
Z can be uniquely determined within measure- 
ment fluctuations. There are two complications 
to these simple relations: the signal output of the 
scintillator does not scale exactly like Z^, and 
contributions of signals due to J-rays lead to an 
additional component in the Cerenkov light and 
generate a slight increase in the scintillation signal 
above the minimum ionization level. These effects 



have b e en discussed previously bv lCahbauer et al" 



(|2004D : lAve et al.l (|2008ll201li . 

The energy measurement of TRACER depends 
on the energy response of the Cerenkov signal, the 
magnitude of the relativistic increase in specific 
ionization in the gas-proportional tubes, and on 
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Fig. 3. — Schematic of the dual gain readout of 
the proportional tubes. 

the energy dependence of the transition radiation 
(TR) signal. The characteristic response curves 
are shown in Figure O The very rapid rise of the 
Cerenkov signal permits an accurate measurement 
at low energies (■^O.S to 3 GeV amu~^). The much 
slower rise of the ionization signal allows for an 
energy determination in the range of '-^ 10 to a 
few hundred GeV amu~^, if the magnitude of this 
rise is larger than the level of fluctuations in the 
signals. Finally, the highest energies can be rather 
precisely determined from the rapid increase of the 
TR intensity. Because of the steeply falling energy 
spectrum of cosmic rays, the measurement at high 
energies is limited by counting statistics and not 
by saturation of the detector response. 

For LDB2, detector upgrades were imple- 
mented in three areas to include the light elements 
boron (Z = 5), carbon (Z = 6), and nitrogen 
{Z = 7) in the measurement: 

1. The signals from each wire of the propor- 
tional tubes were split resistively into two 
amplifiers (see Figure |3]), increasing the over- 
all dynamic range of from 10^ to ~ 10'' 
(i. e. charge signals from ^ 10^^^ Coulomb 
to nearly 10""'^ Coulomb). This change re- 
quired a doubling of the number of ampli- 
fier ch annels (provided by AMPLEX VLSI- 
chips ( Beuville et al. 1990)), and a complete 
redesign of the digitization and data acqui- 
sition electronics. 

2. The xenon content in the gas mixture for the 
proportional tubes was increased by almost 
a factor of 4 as compared to LDBl. As a 
result, the energy resolution is considerably 
enhanced (see Figure [3] and Section 13. 3p . 
The increased xenon content also causes a 
slight change in the TR response. The up- 
grade facilitates the energy measurements 
for boron and carbon in the 100 GeV amu^'^ 
energy range. 




Fig. 4. — Flight trajectory of LDB2 from Kiruna 
(Sweden) to Northern Canada. 



3. While for LDBl only a scintillator on top 
and a scintillator- Cerenkov counter pair on 
the bottom were employed, two identical 
scintillator-Cerenkov counter pairs were now 
incorporated on the top and bottom. This 
permits independent charge measurements 
on top and bottom and results in an overall 
improvement of the charge determination. 

2.2. Balloon Flight 

The TRACER instrument was launched on a 
long-duration balloon flight in the Northern hemi- 
sphere from Kiruna (Sweden) on July 8, 2006. 
While it was intended to complete at least one 
circular path around the North Pole, nominally 
taking 14 days, the flight was terminated over 
Northern Canada after 5 days as permission for 
overflight of Russian territory could not be ob- 
tained. Figure |4] shows the flight trajectory. The 
balloon floated at an average altitude of ^ 38 km 
(corresponding to 4.5 g cm~^ of vertical resid- 
ual atmosphere), with daily variations as shown 
in Figure \5\ The instrument operated well, and 
~ 3 X 10^ cosmic-ray events were recorded. All 
electronics performed as expected. The propor- 
tional tube system performed without problems 
during the whole flight. After 1.5 days of flight, a 
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Fig. 5. — Profiles of altitude and vertical atmo- 
spheric overburden as a function of flight time for 
LDB2. 



Fig. 6. — Two-dimensional charge distribution for 
light nuclei. 



subset of PMTs in the bottom scintillator and in 
the top Cerenkov detector had to be deactivated 
due to an electric discharge. This had only a minor 
impact on the instrument's performance. 

For each event, the signals from all PMTs are 
recorded. A zero-suppression system is utilized 
for the proportional tubes, and only signals above 
a threshold value are stored. All events are as- 
sembled and processed by the on-board CPU for 
storage on a hard disk, and are also sent to the 
telemetry transmitters. A total of 50 GB of raw 
data, at an event rate of ^ 120 Hz, were gathered. 
The lifetime of the instrument was 203,200 s, at a 
dead time of the data acquisition of ^ 18%. 

3. Data Analysis 

The data analysis follows the struct ure of the 
analy sis performed for the data of LDBl (jAve et al 
20081 ). Events are accepted only if the follow- 
ing "quality cuts" are satisfied: (1) a minimum 
of three PMT's in each of the scintillation and 
Cerenkov counters must show a non-zero signal; 
(2) the event must correspond to a well recon- 
structed trajectory; and (3) the signal amplitudes 
measured in the 16 proportional tube layers must 
be consistent within expected levels of fluctua- 
tion. These cuts have a combined efficiency of 
^ 80%, independent of energy (see Table [1]). An 



outline of the analysis is given in the following, 
with emphasis on the effect of the modifications 
for LDB2. 

3.1. Trajectory Reconstruction 

First, the trajectory of every cosmic-ray nucleus 
through the instrument is determined. A first- 
level reconstruction that fits straight lines to the 
center positions of all tubes that recorded a signal 
yields a lateral precision of ^ 5 mm. A second- 
level fit utilizes the fact that the relative signal 
amplitudes are proportional to the track length in 
each tube, and thus are a measure for the "im- 
pact parameters" of the track from the tube cen- 
ter. The track fit then has a lateral precision of 
^ 2 mm. This corresponds to an uncertainty in 
total track length through the proportional tubes 
of - 3%. 

3.2. Charge Measurement 

The charge measurement relies on the combina- 
tion of the signals of the scintillat or and Cerenkov 
detec tors as described previously ( Ave et al. 2008L 
20111 ). The additional Cerenkov detector used 
in LDB2 facilitates independent charge measure- 
ments at both the top and bottom of the instru- 
ment. 
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Fig. 7. — Parametrization of the charge distribu- 
tion of top and bottom charge for hght nuclei. The 
shape of the selection cut is indicated. The on and 
off region method to estimate the carbon contam- 
ination in the boron data sample is illustrated. 



The charges obtained for each event are shown 
in Figure [6] for the light elements boron, carbon, 
nitrogen, and oxygen. The elements are clearly 
separated with an effective charge resolution of 
0.23 charge units for boron and carbon nuclei, and 
0.25 charge units for nitrogen and oxygen nuclei. 
The resolution rises for higher charges, reaching 
0.55 charge units for iron nuclei. 

The distributions in Figure [5] are parametrized 
with two-dimensional Gaussian functions and are 
shown in Figure [71 Semi-elliptical cuts, commen- 
surate with the resolution of the distributions, are 
applied as indicated in the figure. The cut val- 
ues Rtop and Rbot (see Figure [T]) are chosen to 
minimize the contamination by adjacent charges 
and are different for the individual elements. For 
instance, for iron nuclei Rtop — ±1.0 (in charge 
units) while for boron nuclei, Rtop — ±0.3. The 
cut Rbot is asymmetric, with Rtot = (—1-2/ + 2.4) 
for iron, but Rtot = (—0.32/ -I- 0.64) for boron. 
The abundant elements oxygen and iron are then 
selected with an efficiency of 68% (see Table [T]) 
with negligible contamination from neighboring el- 
ements. Tighter cuts applied to boron and carbon 
nuclei lead to efficiencies of 36% and 44%, respec- 




15 20 25 30 

bottom charge Z 

Fig. 8. — Bottom Charge Zi,ot distribution of 
events identified as iron at the top of the in- 
strument. The solid line represents a Gaus- 
sian fit, which illustrates the distribution of non- 
interacting nuclei. 



tively. Even then, the selected boron events at 
high energy contain a very small carbon contami- 
nation of 0.63 ± 0.03%. This contamination level 
has been determined by two independent methods: 
(1) by moving the boron selection cut around the 
carbon peak to an "off" position (see Figure [7]) 
and counting the events within the cut; (2) by in- 
tegrating the tail of the parametrized carbon dis- 
tribution within the selection cut for boron. Both 
methods lead to the same result. There is no con- 
tamination of the boron sample by beryllium nu- 
clei. The effect of the contamination of the boron 
events on the determination of the boron energy 
spectrum is discussed in Section H?^ 

A fraction of events undergo a charge-changing 
interaction within the material of the instrument 
(7.7 g cm~^). Such events exhibit an apparent 
lower nuclear charge at the bottom of the instru- 
ment, since the sum of the signals of the daughter 
nuclei is smaller than the signal of the parent nu- 
cleus. The interacting fraction can either be ex- 
perimentally determined, as illustrated in Figure|8] 
for iron nuclei, or can be calculated, using the 
in teraction cross sections i n the parametrizatio n 
of' Bradt and Peters! (|l950l ): IWestfall et al.l (Il979l) . 
which are assumed to be independent of energy. 
The results of both methods are in agreement and 
one obtains a survival probability for iron nuclei 
through the instrument of ^ 50%. Similarly, the 
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Bradt and Peters formula is used to determine in- 
teraction losses in the residual atmosphere above 
the balloon. The survival probabilities for all ele- 
ments considered in the present work are given in 
Table HJ They are evaluated at an average angle of 



incidence of 6* = 30° (see also Ave et al. ( 201ll )) 



The interaction losses are taken into account in the 
calculation of the effective aperture of TRACER 
for each element (see Table [IJ . 

3.3. Energy Assignment 

The energies, or the Lorentz-factors 7 = 
E/mc^ , of the individual primary cosmic-ray nu- 
clei (C, O, Fe) are determined by the same pro- 
cedures that were e mployed for LDB l and were 
described previously ( Ave et a"l1l2008l ). 

At low energies (^0.5 to 3.0 GeV amu"^), the 
energy measurement is performed from the signals 
of the bottom Cerenkov counter. The counter has 
a steep energy response in this region (see Fig- 
ure El) , hence, it exhibits excellent energy resolu- 
tion. 

At high energies, the dE/DX signals from the 
entire proportional tube array provide the en- 
ergy measurement over the range of ^ 10 to 
400 GeV amu~^. The correlation between the 
dE/dx-array and the TRD identifies particles 
above ~ 1000 GeV amu~^. Because of the sudden 
change in the response function around the onset 
of TR, the energy resolution is not well defined in 
the energy range of about 400 to 1500 GeV amu^^. 
Therefore, this energy region is excluded from 
analysis. 

The increase in the xenon content of the pro- 
portional tubes modifies the detector response in 
two respects: (1) the relative magnitude of the 
relativistic increase in dE/dx is slightly reduced 
to 33% instead of 40% over the 7-range from 10 
to 440. However, the absolute signal is increased 
by nearly a factor of 4. Consequently, the relative 
signal fluctuations are reduced by about a factor 
of two, and the relative energy resolution between 
10 and several 100 GeV amu~^ is improved by the 
same factor. This important fact is illustrated in 
Figure ini (2) for the TRD, both the signal contri- 
bution from ionization, and the conversion proba- 
bility of x-rays are increased as compared to LDBl 
because of the larger xenon content of the tubes. 
In balance, the appearance threshold for x-rays 
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Fig. 9. — Measured energy resolutions {la) for 
LDBl and LDB2 for the dE/dx-array. 



moves to a higher 7-value of 785 ± 140 and the 
relative magnitude of the TR signal decreases by 
19 ± 3%. More deta il on these calibrations is given 
bv lAve et al. (2011) and references therein. 

The rather steep energy dependence of the TR 
response curve (Figure [2]) insures good energy res- 
olution for the highest energy events. The mea- 
sured energy resolution for all three sub-detectors 
is shown in Figure [10] as a function of charge Z. 

The accepted high-energy events are character- 
ized by two quantities, the average ionization loss 
(dE/dx), and the average TRD signal [dE/dx + 
TR), which are defined as the sums of the sig- 
nals in the dE/dx-array and TRD, respectively, 
divided by the sum of the track lengths. It is vital 
for the energy assignment in the high energy re- 
gion that particles at energies below the minimum 
ionization level are excluded as the response func- 
tion (Figure [2]) is double- valued and rises towards 
lower energies. This is accomplished with high ef- 
ficiency by a cut on the bottom Cerenkov signal 
at the level corresponding to minimum ionization 



(7 = 3.97). 

Figure [Til shows a cross correlation of (dE/dx) 
and (dE/dx + TR) for oxygen nuclei (Z=8) above 
3 GeV amu~^. At energies below TR threshold 
(~ 1.4 MIP, where MIP is the signal of a minimum 
ionizing particle, at an energy of 3 GeV amu^^) 
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Fig. 10. — Measured energy resolution {la) of in- 
dividual detector subsystems vs. charge Z for typ- 
ical energies for LDB2. 



both signals are the same within fluctuations. The 
distribution peaks around 1 MIP and rises along 
the diagonal as the energy increases. After the on- 
set of TR, the events rise above the diagonal with 
increasing energy. Unambiguous TR-events, well 
outside of possible dE/dx fluctuations, are high- 
lighted in the scatter plot. These events have ener- 
gies well above 1,000 GeV amu~^. Note that their 
number represents less than 10"'* of the total num- 
ber of events in the plot. Nevertheless, the plot is 
remarkably free of background and demonstrates 
the power of identifying particles by independent 
measurements of ionization loss and TRD signal. 
The dashed lines in the plot also indicate the edges 
of the energy bins into which the data are sorted. 
The same procedures have been followed for the 
primary elements iron (Z=26) and carbon (Z=6). 

To generate the energy spectra of the individual 
elements, the events are sorted into energy bins 
as indicated (with a gap between about 400 to 
1500 GeV amu~^ for reasons of energy resolution), 
and "overlap corrections" are applied to account 
for possible misplacements into neighboring bins. 
The corresponding correction factors are given in 
Table [TJ They are close to unity as the energy 
bins are chosen to be quite wide with respect to 
the energy resolution. The exception are the low 



Fig. 11. — dE/dx-array signal vs. TRD signal for 
high-energy oxygen data. The diagonal and en- 
ergy bin edges (dashed) are indicated, as well as 
the TR rsponse (dash-dotted). Minimum ionizing 
energy is indicated as a white marker. Possible 
TR events are shown in bold symbols. 

energy dE/dx bins for boron and carbon, where 
signal fluctuations are large. 

4. Energy Spectra 

4.1. The Primary Nuclei C, O, and Fe 

The differential intensity for each element in an 
energy bin j is derived from the number of events 
ANj in the bin of width AEj as: 



dN\ _ 1 AN J ujj 



dE 



(1) 



with the total exposure $7, the total efficiency 
Cj, and the overlap correction ujj. For a power- 
law spectrum with spectral index a, the in- 
tensity is plotted a,t an energy Ej, defined as 
(|Laffertv and WvattlflQQSh : 



1 



AEi 1 



E. 



E' 



-I/O 



(2) 



with AEj — Ej+i — Ej. Spectral indices a com- 
mensurate with previous measurements are used 
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Table 1 

Summary of efficiencies and exposure factors. 



Item 


Spectrum Range 


All Elements 


Boron 


Carbon 


Oxygen 


Iron 


survival probability 


all 




U.I 6 


n 70 


U.Dy 


U.4o 


instrument 














survival probability 


oil 

an 




U.oU 


n 7Q 


U. ( / 


n fit; 


atmosphere 














exposure factor 


oil 

aii 




O.ZZ 


K no 


A VO 
4. / z 


Q1 

z.oi 


(m^sr days) 
















CER 


0.043 










zenith cut 


dE/dx 














TRD 














CER 




0.78 


0.79 


0.79 


0.80 


quality cuts 


dE/dx 




0.79 


0.80 


0.81 


0.79 




TRD 




0.79 


0.80 


0.81 


0.79 




CER 




0.49 


0.56 


0.81 


0.78 


charge cut 


dE/dx 




0.36 


0.44 


0.68 


0.67 




TRD 




0.36 


0.44 


0.68 


0.67 




CER 




0.9 


0.9 


0.9 


1.0 


correlation norm. 


dE/dx"- 




0.2/0.7 


0.6/0.9 


0.7/0.9 


0.9 




TRD 




1.0 


1.0 


1.0 


1.0 



"Two individual values are given for the two dE/dx energy bins since they differ significantly. 



(2.65 at high energies for primary elements). How- 
ever, the value of Ej is not very sensitive to the 
choice of a. The highest-energy data point is an 
integral value above Ej, i. e. — ^ c». The 
corresponding differential intensity is then 



dN 
'dE 



AN, 



1 

n 



Wo 



(a - 1)E 



— a 
■med 



E 



and is plotted at the median energy Emed' 



(3) 



(4) 



The differential energy spectra for iron, oxygen, 
and carbon nuclei are given in terms of absolute 
intensities on top of the atmosphere in Figure [T2l 
and are summarized in Table [2] Uncertainties 
are stated separately for statistical and systematic 
contributions. The figures include previously pub- 
lished data for comparison, and the energy spectra 
are shown multiplied by _e2-65 j-qj. d^j-ii^y^ 

The statistical uncertainties are based on Pois- 
son counting statistics. For ten or fewer events, er- 
rors become asymmetric and take values according 



to iGehreld (|l986t ) for a confidence level of 849^0. 
The systematic errors represent uncertainties for 
energy bin overlaps, selection cut efficiencies, and 
imperfect knowledge of the response functions. 

4.2. Energy Spectrum of Boron 

Boron is a secondary element produced by spal- 
lation mainly of carbon and oxygen. These parent 
nuclei are much more abundant than boron nuclei. 
Therefore, two additional considerations must be 
taken into account when deriving the energy spec- 
trum of boron nuclei: (1) contamination of the 
boron sample with a small admixture of carboir 
nuclei (see Section !??^ . and (2) the contribution of 
boron nuclei generated in the residual atmosphere 
above the balloon. 

To investigate the carbon contamination, we 
consider the scatter plot {dE/dx + TR) vs. 
(dE/dx) shown in Figure [T31 This plot is derived 
from about 40,000 boron nuclei. It shows qualita- 
tively the same behavior as Figure [TT] for oxygen. 



^More recent work bv lFeldman and Cousin"i l ll99St) leads to 
very similar values. 
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Table 2 
Measured intensities. 



Energy Range 


Energy E Number of Events 


Intensity ±(Tstat ± fsys 


GeV amu~^ 


GeV amu~^ 


(m^ sr s GeV amu~^)~^ 



Iron (Z=26) 












0.8 - 1.0 


0.9 


593 


(4.1 


± 0.2 ± 0.1) 


X 10"^ 


1.0 - 1.3 


1.1 


487 


(2.78 


± 0.1 ± 0.09) 


X 10"^ 


1.3 - 1.7 


1.5 


386 


(1.79 


± 0.07 ± 0.2) 


X 10"^ 


14-38 


23 


1180 


(4.8 


± 0.1 ± 1.3) 


X 10"* 


38 - 101 


60 


314 


(4.7 


± 0.2 ± 1.7) 


X 10"^ 


101 - 264 


159 


57 


(3.2 


± 0.3 ± 1.4) 


X 10"^ 


1200 - 


1800 


1 


(5 


I 10 ,2\ 
±4 ±3) 


X 10"^ 



Oxygen (Z=8) 



0.8 - 1.0 


0.9 


7900 


(3.26 


± 


0.04 ± 0.09) X 


10" 


1.0 - 1.3 


1.1 


6440 


(2.17 


± 


0.03 ± 0.05) X 


10° 


1.3 - 1.7 


1.5 


6210 


(1.28 


± 


0.02 ± 0.05) X 


10° 


9.5 - 87 


26 


43700 


(2.25 


± 


0.01 ± 0.06) X 


10"^ 


87 - 432 


181 


2140 


(1.91 


± 


0.04 ±g:i) X 


10"'^ 


1500 - 


2300 


12 


(1.9 


± 


0.6 ±?f) X 


io-« 



Carbon (Z=6) 



0.8 - 1.0 


0.9 


6510 


(3.19 ± 0.04 ± 0.1) X 10° 


1.0 - 1.3 


1.1 


5460 


(2.19 ± 0.03 ± 0.1) X 10° 


1.3 - 1.7 


1.5 


4760 


(1.19 ± 0.02 ± 0.1) X 10° 


9.5 - 87 


26 


28300 


(2.17 ± 0.01 ± 0.06) X 10" 


87 - 432 


181 


1430 


(1.70 ± 0.04 ± S I) X 10" 


1500 - 


2300 


10 


(2.2 ±°:?±°:i) xio- 



Boron^ (Z=5) 



0.8 - 1.0 


0.9 


2170 


(1.13 


± 0.03 ± 0.1) 


X 


10° 


1.0 - 1.3 


1.1 


1600 


(7.3 


± 0.2 ± 0.1) 


X 


10" 


1.3 - 1.7 


1.5 


1280 


(3.8 


± 0.1 ± 0.1) 


X 


10" 


9.5 - 86.5 


23.8 


7200 (27) 


(4.9 


± 0.1 ± 1.0) 


X 


10" 


86.5 - 432 


173 


413 (48) 


(1.6 


± 0.1 ± 0.4) 


X 


10-' 


1500 - 


2070 


5(4) 


(4 


±3 ±2) 


X 


10- 



*in brackets: expected number of carbon background 
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Fig. 12. — Differential energy spectra of car- 
bon, oxygen, and iron (top to bottom). The 
spectra are multiplied by E^-^^, as well as di- 
vided by 10 for oxygen and iron. For cqmpar i- 
son, re sults from HEAO (Engelmann et al. 1990f). 



CRN (iMiiller et all 19911). AT IC 



Panovetal 



20071 ). and CREAM (jAhn et al.ll2009l) are shown. 
The error bars are statistical only. 



although the relative statistical fluctuations are 
larger due to the lower Z of boron. However, an 
analysis of the charge selection for boron (see Sec- 
tion 13. 2|) based on Figures |6] and [71 predicts that 
FigurelTSlmust include a contamination of 250±12 
carbon nuclei (0.63 ± 0.03%), clustered around 
the MIP signal of carbon events. The effect of 
this background is illustrated in the histogram of 
{dE /dx+TR) for boron in Figure [Ml which repre- 
sents the projection of the scatter plot of FigurefTSl 
upon the vertical axis. As expected, the histogram 
is asymmetric towards higher signals due to the 
relativistic rise and due to the appearance of TR 
photons. Also included in Figure [14] is a simulated 
distribution of the 250 carbon background events. 
This background resides mostly in the region of 
tail of the boron distribution, peaking near the 
onset of TR. Obviously, the carbon contamination 
does not affect the boron measurement at low 
energy, below ^ 100 GeV amu~^. However, the 
contamination makes the identification of boron 
unreliable in the region 450 — 1500 GeV amu~^. 



N 

a: 



22.5 

in 
Q 

E 2 



1.5 




1 1.5 2 2.5 3 

dE/dX signal [MIP/Z^] 

Fig. 13. — dE/dx-array signal vs. TRD signal for 
high energy boron data. Small and thick black 
markers represent events at low and high energy, 
respectively. Minimum ionizing energy is indi- 
cated as a white marker. The diagonal is shown as 
a dashed line, and the dash-dotted line represents 
the TR response. 



which is already excluded from analysis due to 
the ill-defined energy resolution around the onset 
of TR. There still remains at least one event at 
very high energy which very unlikely (at a prob- 
ability of 10~^) is a carbon nucleus. This event 
also clearly stands out in the scatter plot of Fig- 
ure [13] and appears to be a genuine boron nucleus 
of ~ 6, 000 GeVamu-i. 

The second effect to be accounted for is the gen- 
eration of secondary boron nuclei by spallation of 
heavier nuclei in the residual atmosphere above 
the balloon. Each galactic component is atten- 
uated in the atmosphere, at an average column 
densitjH X = 5.2 g cm~^, by a factor exp(— X/A). 
Here, A = to/ct is the attenuation path length, 
with m the average mass of the target atoms in 
air, and a the spallation cross section. This atten- 
uation (see also Section 13.21) has been taken into 
account in deriving the energy spectra of the pri- 
mary elements (Section |4?T|) . 



^Averaged over the duration of the flight and over the zenith- 
angle distribution of accepted particles. 
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10 



TRD signal [MIP/Z"] 



Fig. 14. — Histogram of TRD signals. The solid 
histogram represents boron data. The thin line in- 
dicates the distribution of minimum ionizing nu- 
clei. The dashed histogram shows a distribution 
of contaminating carbon events as expected from 
simulation. Note that the vertical axis is logarith- 
mic. 

The contribution of spallation products to the 
observed particle flux must be included for rare 
elements such as boron: the measured intensity 
Nb{X) is then 

Nb{X) = 7Vi3(0)-exp(-X/AB)+^7Vp(0)X/Ap_s 

(5) 

Here, Np refers to parent nuclei heavier than 
boron, and Ap_B is the differential spallation path 
length in air for an interaction of the parent lead- 
ing to boron nuclei. To obtain Nb{0), the boron 
intensity at the top of the atmosphere, this ex- 
pression has been evaluated using differential cross 



sections from Webber et al. I (ll990l) . Contributions 



to the sum in Eq. ([5]) come essentially only from 
carbon and oxygen, and, based on the results de- 
scribed in Section l473l the carbon-to-oxygen abun- 
dance ratio is taken to be independent of energy. 

An independent technique to determine the at- 
mospheric boron contribution utilizes the daily 
variations of atmospheric depth during the flight 
(see Figure [5]) . The corresponding variations in 
the counting rates for the individual elements 



3 

E 
C5 

CO 



3 
E 

IS 

a 



10- 
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10 



^ I I J 
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Fig. 15. — Differential energy spectrum of boron. 
The spectrum is multiplie d by E^-^^ . For compar - 
ison, result s from HEAP jEngelmann et ah 19901 ) 
and CRN (jSwordv et all Il990l ) are shown. The 
error bars are statistical only. 



are a direct measure of the atmospheric effects. 
The analysis leads to results which are consistent 
with those obtained with published cross sections. 
More detail on this procedure will be given in a 
separate publication dOb ermeier et al.l 120111) . 

The contribution of atmospheric secondary 
boron nuclei to the measured energy spectrum 
is small at low energy (~ 6% at 1 GeV amu~^), 
but reaches ~ 20% above 1 TeV amu~^. Thus, 
the single boron event defining the highest-energy 
data point has a 20% chance of being of atmo- 
spheric origin. The energy spectrum of boron af- 
ter corrections for atmospheric contributions and 
interaction losses in both, atmosphere and instru- 
ment, is shown in Figure [TSl and Table [21 



4.3. Relative Intensities of Primary and 
Secondary Cosmic Rays 

It was observed in LDBl that the primary spec- 
tra of the primary cosmic-ray nuclei from oxygen 
to iron above 20 GeV amu~^ can be described by 
a power law in energy with a common spectral 
index of 2.67 ± 0.05. Consequently, the relative 
abundances of the different elemental species do 
not change with energy. The same behavior char- 
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Fig. 16. — The carbon-to-oxygen abundance 
ratio as a function of kinetic energy per nu- 
cleon. Error bars are statistical only. Prev i- 
ous res ults from HEAP dEngelmann et aP 19901) 



CRN (ISwordv et al.l Iwm. ATIC (IPanov et al 



2007bl ). and CREAM (jAhn et al.ll2008l) are shown 
for comparison. 

acterizes the results of LDB2, including the pri- 
mary element carbon. This is iUustrated in Fig- 
ure [16] which shows the abundance ratio carbon to 
oxygen as a function of energy. This figure also 
includes previous results from other experiments. 

The relative abundance of boron, which is of 
prime importance for the understanding of Galac- 
tic propagation of cosmic rays, is expressed by the 
boron-to-carbon {B/C) abundance ratio. This ra- 
tio is calculated from the differential intensities in 
Table [21 and presented in Table [3| and Figure [iTj 
as a function of kinetic energy per nucleon. While 
the same energy intervals have been used to de- 
rive the energy spectra, the mean energy in each 
bin differs slightly because of the different spectral 
slopes (Equations ([2]) and (|4|)). For calculating the 
ratio, the carbon flux has been scaled to the same 
energies as are used for the boron spectrum. The 
table states both statistical and systematic uncer- 
tainties. In the figure, statistical and systematical 
uncertainties are shown as thin and thick error 
bars, respectively. 

The statistical uncertainty of the boron-to- 
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CRN 
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Level of atmospheric B production 
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Fig. 17. — The boron-to-carbon abundance ra- 
tio as a function of kinetic energy per nucleon. 
Error bars are statistical (thin) and system- 
atic (thick). A model corresponding to an es- 
cape path length oc E~^-^ (dotted) is shown. 
The level of atmospheric production of boron 
(dashed) is indicated. Previous measurement s 
are sh own from HEAP (Engelmann et al. 1990[) 



CRN (ISwordv et al 



1990^ 



ATIC (Panov et al 



2007bl). CREAM lAhn et all |2008[) and AMS 



01 ([Aguilar et al.ll2010[) . 



carbon ratio is de rived with a pr ocedure based on 
Bayes' Theorem ( Paternol 2003 ). which extends 
Poisson statistics to small numbers and small ra- 
tios. The dashed line in Figure [ITj indicates the 
expected level of the boron-to-carbon ratio due 
to atmospheric production of boron nuclei, which 
has been subtracted. 

Many systematic uncertainties, i.e. for aper- 
ture, efficiencies, and detector response functions, 
are correlated for the measurement boron and car- 
bon nuclei, and thus cancel for the ratio. The main 
source of systematic uncertainty remaining comes 
from the calculation of the amount of boron nu- 
clei produced in the atmosphere. Also, the uncer- 
tainty of the energy measurement does not entirely 
cancel due to the different spectral indices of the 
boron and carbon energy spectra. 
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Table 3 
The boron-to-carbon ratio. 



Energy Range Energy E 
GeV amu~^ GeV amu~ 



0.8 - 1.0 0.9 

1.0-1.3 1.1 

1.3-1.7 1.5 

9.5 - 86.5 23.8 

86.5 - 432 173 

1500 - 2070 



B/C ratio ±astat ± ""s; 



( 3.6 ± 0.1 ±8:1 ) X 10"^ 

( 3.3 ± 0.1 ±EI:i ) X 10"^ 
( 3.2 ± 0.1 ±8:i ) X 10'^ 
( 1.80 ± 0.04 ±El:[3i ) X 10^ 

( 8.2 ± 0.9 ±li ) X 10-2 

( 1.2 ±Ji±8:t ) X 10-1 



5. Discussion and Conclusions 

We have described a new measurement of the 
energy spectra of highly relativistic cosmic-ray nu- 
clei. The TRACER instrument had previously 
(LDBl) provided measurements of the heavier nu- 
clei from oxygen to iron (Z = 8 to 26). The most 
recent flight (LDB2 in 2006) includes the lighter 
nuclei and covers the region from boron to iron 
(Z = 5 to 26). LDBl had a longer duration, and 
therefore better statistics than LDB2. Neverthe- 
less, the two measurements show absolute intensi- 
ties that agree with each other over the wide range 
of energy covered. This provides confidence that 
the new results for the elements carbon and boron 
(Figures [14] and [T5\\ are consistent with the previ- 
ous work. 

It should be noted, however, that the differ- 
ential intensities from the present work at low 
energy (at about 1 GeV amu"^) are consistently 
higher compared to previous measurements. This 
can be attributed to the differe nt solar modula- 
tion of low-energy cosm i c ray s (jGoldstein Il970t 
Buchvarova and Velinov 2010f) at the time of 
the flight. Where overlap exists, the TRACER 
results are also in goo d agreement with mea - 
sureme nts from HEAP jEngelmann" et al.l 1X99(^1 



CRN (Mulleretal 



1991 



ATIC ( Panovetal 



2007t ). and CREAM lAhn et al.| [2009t). but have 



better statistics than some of these measurements. 
Currently, the analysis of the TRACER LDB2 
data is still ongoing, concentrating on the spectra 
of the primary elements neon (Z = 10), magne- 
sium {Z = 12), and silicon {Z = 14). Results on 
these elements will be published in due course. 
A compilation of all present results from the two 
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Fig. 18. — Compilation of the differential en- 
ergy spectr a measured by T RACER in LDBl and 
LDB2 (see lAve et all (|2008l ) and this work). The 
dashed lines represent a simple power-law fit above 
20 GeV amu-i. 



LDB flights of TRACER is given in Figure [HI 

For boron, there are only two measurements 
besides TRACER that report a bsolute intensities 
at re lativistic en ergies: HEAP (jEngelmann et aT 



1990h and CRN (ISwordv et al.lll990h . but none of 
these extends into the TeV amu"^ region. These 
data are included in Figure [151 Good agreement 
with the present work can be noted within the 
quoted error limits, and in the energy region where 
overlap exists. 

Most important is the energy dependence of 
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the relative boron intensity, i. e. of the B/C ra- 
tio. In Figure [171 the TRACER results are com- 
pared t o previous work, inc ludinR the r esults from 
ATIC (jPanovet al. l l2007bl). CREAM dAhn et al ' 



20081) . and AMS-01 (lAeuilar et all 120101) . for 



which the absolute energy spectrum of boron has 
not been reported. Figure [17] also includes an 
extrapolation (dotted line) for an E~^-^ decrease 
of this ratio with energy o n the basis of mea- 
surements at lower energv (lYanasak et aL 2001 



Engelmann et al.lll990t Swordv et al .I ll990l) . How 
ever, there is no model of propagation in the ISM 
that would predict such a behavior a priori. Fur- 
thermore, this energy dependence would imply un- 
reasonably small values for the propagation path 
len gth of cosrn i c rays at the highest energies (see, 
e.g. [Ave et aL I (l2009l) l. Thus, it is quite interest- 
ing that the data points of TRACER on the B/C 
ratio above 100 GeV amu~^, and also, as shown 
in Figure [T71 some of the results of other measure- 
ments, appear to lie above the E~^-^ prediction. 
This feature may have profound implications upon 
our understanding of cosmic-ray propagation. It 
may suggest that the energy dependence of the 
propagation path length flattens at high energy, 
and perhaps indicate an asymptotic transition to 
a constant residual path length. A more detailed 
analysis of this situa tion will be presented in a 
separate publication ( Obermeier et al. 2011b[ ). 

In conclusion, TRACER has proven to be an ex- 
cellent detector for measuring the B/C ratio in the 
TeV amu~^ region. The instrument has outstand- 
ing discrimination power to identify boron nuclei 
at the highest energies, i.e. above ~ 1 TeV amu~^, 
where contamination with carbon nuclei is no 
longer a concern. A 30-day flight could yield ~ 10 
boron events above 1 TeV amu~^ and thus, pro- 
vide a definitive measurement of the propagation 
path length at those energies. 

We are grateful to the services of the technical 
staff at the University of Chicago, in particular 
the contributions of G. Kelderhouse, R. Northrop, 
D. Plitt, C. Smith, and P. Waltz. We thank Dr. 
M. Ichimura for help with the data analysis, and 
we acknowledge the dedicated work of a group of 
undergraduate students led by E. Brannon, and 
supported by the Illinois Space Grant. We are in- 
debted to the staff of the Columbia Scientific Bal- 
loon Facility and the Esrange Space Center (Swe- 



den) for the successful balloon flight operations. 
AO acknowledges support of FOM in the Nether- 
lands ("Stichting voor Fundamenteel Onderzoek 
der Materie" ) . This work was supported by NASA 
through grants NNG 04WC08G, NNG 06WC05G, 
and NNX 08AC41G. 

REFERENCES 

Aguilar M. et al., 2010, ApJ, 724, 329 

Ahn H. S. et al., 2007, NIM A, 579, 1034 

Ahn H. S. et al., 2008, Astro. Ph., 30, 133 

Ahn H. S. et al., 2009, ApJ, 707, 593 

Ahn H. S. et al., 2010, ApJ, 714, L89 

Ave M. et al., 2008, ApJ, 678, 262 

Ave M. et al, 2009, ApJ, 697, 106 

Ave M. et al., 2011, NIM A, in print, DOI: 
10.1016/j.nima.2011.05.050 

BeuviUe E. et al, 1990, NIM A, 288, 157 

Bradt H. L. and Peters B., 1950, Phys. Rev., 77, 

54 

Buchvarova M. and Velinov P. I. Y., 2010, Adv. 
Space Res., 45, 1026 

Engelmann J. J. et al., 1990, A&A, 233, 96 

Feldman G. J. and Cousins R. D., 1998, 
Phys. Rev. D, 57, 3873 

Gahbauer F. et al., 2004, ApJ, 607, 333 

Ganel O. et al., 2005, NIM A, 552, 409 

Garcia-Munoz M., Mason G. M., and Simpson J. 
A., 1977, ApJ, 217, 859 

Gehrels N., 1986, ApJ, 303, 336 

Goldstein M. L. et al., 1970, Phys. Rev. Lett., 25, 
832 

Juhusson E., Meyer P., and Miiller D., 1972, 
Phys. Rev. Lett., 29, 445 

Lafferty G. D. and Wyatt T. R., 1995, NIM A, 
355, 541 

L'Heureux J. et al., 1990, NIM A, 295, 246 



15 



MiiUer D. et al., 1991, ApJ, 374, 356 

Obermeier A. et al., 2011, ApJ, in preparation 

Obermeier A. et al., 2011b, ApJ, in preparation 

Panov A. D. et al., 2007, Bull. Russ. Acad. Sci. 
Phys, 71, 494 

Panov A. D. et al., 2007b, Proceedings of the 30th 
International Cosmic Ray Conference (Merida), 
2, 3 

Paterno M., 2003, D0 note 2861 
Swordy S. P. et al., 1990, ApJ, 349, 625 
Webber W. R. et al, 1990, Phys. Rev. C, 41, 533 
Westfall G. D. et al, 1979, Phys. Rev. C, 19, 1309 
Yanasak N. E. et al., 2001, ApJ, 563, 768 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 



16 



